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Abstract —We consider balanced flows in a natural gas trans¬ 
mission network and discuss computationally hard problems 
such as establishing if solution of the underlying nonlinear gas 
flow equations exists, if it is unique, and finding the solution. 
Particular topologies, e.g. trees, are known to be easy to solve 
based on a variational description of the gas flow equations, but 
these approaches do not generalize. In this paper, we show that 
the gas flow problem can be solved efficiently using the tools of 
monotone operator theory, provided that we look for solution 
within certain monotonicity domains. We characterize a family 
of monotonicity domains, described in terms of Linear Matrix 
Inequalities (LMI) in the state variables, each containing at 
most one solution. We also develop an efficient algorithm to 
choose a particular monotonicity domain, for which the LMI 
based condition simplifies to a bound on the flows. Performance 
of the technique is illustrated on exemplary gas networks. 

I. Introduction 

Natural gas is the fastest growing component of the energy 
industry. In the USA this growth is primarily due to the 
hydro-fracking revolution [?], which has provided tremen¬ 
dous increase in supply, driving down prices and stimulating 
significant network expansion and creating new uses of the 
natural gas [?], [?]. The revolution also makes operating, 
controlling and optimizing the emerging gas systems more 
challenging. Local considerations which prevailed in the past 
are no longer sufficient to keep operations of the natural 
gas system economic and reliable. Developing new methods 
to analyze the state of the growing network efficiently and 
accurately has now become a top priority for gas system 
operators and other entities, such as power system operators 
with many gas fired turbine on balance, interested to monitor 
the status of the consumption/production, network flows 
and pressure profile along the long haul gas transmission 
networks. 

Resolving many of the emerging computational challenges 
in operations and planning of the natural gas networks 
depends on solving the so called Gas Flow (GF) equa¬ 
tions. These equations describe spatio-temporal distribution 
of the mass flows and pressures in pipes over an arbitrary 
network subject to known profile of external parameters 
characterizing injection/consumption and compression, see 
e.g. [?], [?], [?]. The equations are nonlinear and difficult 
to analyze, even in the stationary (time-independent) regime 
discussed in this manuscript, assuming overall balance of 
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production and consumption, where thus no accumulation or 
loss of gas in the system (no, so-called, line pack) takes 
place. Analyzing this system of the GF equations entails 
to establishing solution existence, finding the solution and 
establishing uniqueness, if possible. We are not aware of 
publications rigorously analyzing existence and uniqueness 
of solutions to the GF equations. Thus, a natural folklore 
conjecture would be that in its full generality the problems 
of existence, uniqueness and finding a solution are difficult, 
i.e. the problems are of the complexity exponential in the 
size of the network. On the other hand, in a couple of 
special cases the problems were characterized as easy (of 
the polynomial complexity). Special cases are known where 
the GF equations can be solved easily: For tree networks, the 
GF equations can be solved by dynamic programming [?], 
[?], [?]. Another easy case is the one of an arbitrary loopy 
network with no compression, or including only additive 
(and not multiplicative) compression, when the problem of 
finding the balanced gas flow solution becomes equivalent to 
solving a convex optimization problem, that of minimizing 
cumulative losses in the gas pipes due to turbulent friction 
subject to flow conservation at any junction of the network 

m, [?]. 

In this manuscript, we apply the tools of monotone 
operator theory to analyze the GF equations. Just like a 
convex optimization problem is solvable in polynomial time, 
zeros of monotone operators can be found in polynomial 
time [?]. Thus, if we can show that the GF equations can 
be written in terms of a monotone operator, they can be 
solved efficiently. In this manuscript, we indeed show that GF 
operator is monotone over restricted monotonicity domains. 
Each monotonicity domain can be expressed as a Linear 
Matric Inequality (LMI) in the system state (the pressures 
and flows in the network). The monotonicity domains are 
parameterized by a matrix-valued parameter. Within each 
monotonicity domain, there can be at most one solution of 
the GF equations and determining existence and computation 
of solutions can be done in polynomial time. We then show 
that for a particular choice of the monotonicity domain, 
the LMI condition can be replaced with a simple bound 
on the ratio of the maximum to minimum flow in the 
system. Furthermore, the parameters of this monotonicity 
domain can be computed efficiently by solving a quasi- 
convex optimization problem. 

The rest of this paper is organized as follows. In Section 
PH we introduce the GF equations and formulate them as 
a system of multivariate quadratic equations in the flow and 
potential variables. We also introduce the theory of monotone 


operators and solution of variational inequalities. In Section 


in Section [V] we summarize our contributions and discuss 
directions for future work. 

II. Background and Notation 

A. Notation 

R is the set of real numbers. R" denote the corresponding 
Euclidean space in n dimensions. 

Given a set C C R”, Int ( C ) denotes the interior of the 
set. 11 a: 11 refers to the Euclidean norm of a vector x £ R" 
and (x,y) to the standard Euclidean dot product, di {x) is 
the diagonal matrix with with diagonal entries given by x. 
For x, y £ R ra , x * y denotes the vector z with zi = x^. 

Given a differentiable function / : R 1 £ R J , V/ denotes 
the Jacobian of /, a k x k matrix with the *-th row being 
the gradient of the ?'-th component of /. 

For M £ R” xn , Sy(M) = M + M T . S k denotes the 
space of k x k symmetric matrices. 

B. Modeling Gas Networks 

The network is specified by a set of buses (nodes) V and 
a set of gas pipelines (edges). Nodes are denoted by i £ V = 
{0,1,..., n} and edges by (i, j) £ 8. Edges are directed and 
for any pair of nodes only one edge is present, (*, j) £ 8. 
The direction of edges has no physical significance (the flow 
of gas can be directed either way) and it is introduced for 
notational convenience. In some cases we will also denote 
edges by a single index k. Then, k\ is the “head” of the edge 
and &2 is the tail. Let |£| = m. 

In the steady-state, the gas network is characterized by 
pressures at every node and the gas flows over each edge. 
We denote by faj the flow from i to j over edge (i,j) £ 8. 
cf>ij is real and it can be positive or negative. We write i -A j 
if there (i,j) £ 8 and j —> i if (j, i) £ 8. We also denote 
by 7 u the squared pressure at node i. 

The GF equations describe steady-state distribution of 
pressures and flows over the gas network. GF equations 
assume that the network is balanced, i.e. the total sum 
of injections and consumptions is zero. To guaranteer the 
balance one introduces a special node (usually big producer 
or consumer of gas, called the “slack node”) which is flexible 
in its production/consumption. The slack node is denoted by 
0. The remaining nodes are labelled 1 ,,n. Under these 
assumptions the balanced GF equations without compressors 
are given by [?], [?], [?]: 


Kj — Xij (fiij | (pij |, ( 2 , j ) G £ 

(la) 

'y ^ fiij ~ y ] fijii * c {!)■••) 

(lb) 

i-*j j-*i 


7T > 0 

(lc) 


where A,;,- is a known characteristic of the pipe (i,j) 
(constructed from the pipe length, diameter, friction factor 
and the media/gas sound velocity). 


Ill we present our main technical results on the family of 


monotonicity domains of the gas flow equations. In Section 


IV we evaluate our approach on some test networks. Finally, 


Define the nxm matrices A. B. C with entries given by: 


Aik 


Cik = 


1 if k\ = i 

1 if &2 — t 3 ^ik 

0 otherwise 

— 1 if Al2 = z 

0 otherwise 


1 if k\ = i 

0 otherwise 


Then, Eq. ( flb] > can be written as A& = q. 

To solve the GF Eqs. ( |1 a| lb| > means finding the squared 
pressures n, and gas flows, </>, along the pipes of the network 
provided globally balanced injections/consumptions at all 
the nodes and pressure at a node (typically slack bus) are 
known. Note that the PF Eqs. (la B may have no physical 
solution if the pressure set at the slack node is too low, 
i.e. when at least one 7with i > 0 drops below 0. To 
provide an additional pressure boost compressors are intro¬ 
duced. Depending on the local control scheme implemented, 
compressors may be of a multiplicative or additive type. 
Multiplicative compressor boosts the pressure locally by 
constant multiplicative factor in the direction of the flow 
and additive compressor (less standard) provides an additive 
boost along the direction of the flow. 

We assume that the compressors are placed at some 
of the edges in the network, to boost the pressure and 
improve throughput. Compressors are directional. We assume 
(without loss of generality) that the orientation of the edge 
coincides with the compressor direction. 

In the presence of multiplicative compression, the GF 
Eqs. (1 a| lb > generalize to 


Acj) = q (2a) 

O-ij (V., Vjj X jj (pij | |) = (2b) 

7Tj -f“ (1 '('j/j ) A,; .y ©„ | (f>ij V (l ■. j ) £ 8 (2c) 

7T > 0 (2d) 


where r l:l is the relative position of the compressor along the 
edge (i,j) and <*,;,■ > 1 is the compression ratio. We assume 
here that the compressor boosts pressure in the direction from 
i to j. 

It is convenient, for the purpose of further analysis, to 
restate Eqs. (|II-B[) as follows 


A(j> = q (3a) 

(%ij {tti t'ij Ajj(j)ijtpij ) = TTj T (1 Tij ) (3b) 

i’ij = <!>% (3c) 

ipij > 0,7r > 0 (3d) 


where is a newly introduced auxiliary variable. To see 
that Eqs. 0 are equivalent to ( |II-B| ) one observes that 
(|3d|>,(|3c| imply that = \4>ij\- Eqs. (II-Bi constitute a 
set of n + 2 m equations in n + 2?tt, variables: 7r £ R n , </> £ 
£ R m . We order the variables as £ = ( 7 r, t/>), and 

then use Q, to refer to the corresponding blocks of (. 

Eqs. (II-B[) motivates the following definition. 





















Definition 1 . The gas flow operator F : R n+2m R n+2m III. CHARACTERIZATION OF MONOTONICITY DOMAINS 
is defined as of the Gas Flow Operator 


/ A(j>-q \ 

F (£) = I A a 7T — b* (j)*t/} I (4a) 

V J 

where A a = (di (a) B — C) T , b = (a * r + (1 — r))*A > 0. 

Then the GF ( |II-B| ) become simply 

F( C)=0 


C. Monotone Operators 

We now review briefly the theory of monotone operators, 
as it is relevant to the approach developed in the manuscript. 
For details and proofs of the results quoted here, we refer 
the reader to the recent survey [?]. A function H : R fc i —> R fc 
is said to be a monotone operator over a domain C if 

{H(x)-H(y),x-y)>0 Vs.j/eC 

A monotone operator is a generalization of a monotonically 
increasing function (indeed, if k = 1, the above condition 
is equivalent to monotone increase: x > y ==> H (x) > 
H (y)). II is called strictly monotone if 

(H (x) — H (y),x — y) > 0 Vx,y £C,x f=y. 

A common example of a monotone operator is the gradient 
of a differentiable convex function. 

Definition 2 (Monotone Variational Inequality). Let C C R fc 
be a convex set and H be a monotone operator over C. The 
variational inequality (VI) problem associated with H,C is 
given by 

Find x € C such that (H (x), y — x) > 0 My £ C (5) 

We quote the following results from the literature on 
variational inequalities relevant for what follows: 

Theorem II.1. If H is a monotone operator over a convex 
compact domain C and 0- Then, the solution set of the 
variational inequality X* is convex and non-empty. Further, 
an approximate solution x e satisfying 


As discussed in the earlier Sections, the zeros of the 
GF operator correspond to solutions to the GF equations. 
Furthermore, zeros of monotone operators can be found 
efficiently. Thus, if we can prove that the GF operator is 
monotone, then, the GF equations can be solved efficiently. It 
turns out that the GF operator may not be globally monotone, 
but is monotone over restricted domains. We start by defining 
a monotonicity domain formally: 

Definition 3 (Monotonicity Domain). A convex set S C 
R n+2m is said to be a monotonicity domain of the gas flow 
operator F if there exist invertible matrices IV such that the 
operator Fw (x) = WF (x) is monotone over the domain 
x £ S and 

(A 

7T, ip > \phi\,A(f> = q V I cj) £ S 

w 

The motivation behind this definition is that, as long as 
we are interested only in finding zeros of F, it does not 
matter where F or Fw is monotone, since Fw ( x ) = 0 <=$■ 
F (x) = 0 so that zeros of F and Fw are identical. The 
second condition simply ensures that the constraints tt, > 
0, Af = q is satisfied at all points in S, since we require 
these constraints as part of the gas flow equations. 

Theorem III.l. For any invertible matrix W £ 
jgi(n+ 2 m)x(n+ 2 m), convex se t defined by the following 
constraints is a monotonicity domain of F: 


tt > 0,-0 > \fi\,A(j) = q (7) 

/ / ° A 0 \\ 

Sy IV -(di^))- 1 ^ di(V>) di(0) A 0 (8) 

\ V 0 -di(</>) di (VO// 

We denote this set by C (IV). If there is a solution to the GF 
equations in C (IV), it must be a solution to the following 
monotone variational inequality: 


k£C(W), ( F w (k), x - k) > 0 Mx&C(W) 


||x* — x £ || < e for some x* £ X* 


(6) Proof See appendix section [VT| 


□ 


can be found using at most O(d) evaluations of H and 
projections onto C. 

Remark 1. In this manuscript, we will be interested in finding 
zeros of operators corresponding to the GF equations with 
multiplicative compression. It is easy to see that if there 
exists a point x* £ C with H (x*) = 0, then this is a solution 
of the variational inequality. Conversely, if all solutions of 
the variational inequality are such that H (x*) f 0, then 
there is no solution of H (x) = 0 with x £ C. 

Theorem II.2. Suppose H is differentiable. Then H is 
monotone over C if and only if 

XH (x) + XH (xf A0 Mx£C 


Remark 2. The VI stated here may not have a solution since 
C (IV) is not compact. However, for practical gas networks, 
there are always bounds on pressures ir and flow magnitudes 
|<^|. One can add these bounds on tt, 0. y: to the definition 
of C (IV) to get a compact set before solving the VI. Also, 
even when the GF equations have a solution in C (IV), it 
is possible that the VI has multiple solutions (including 
the GF equation solution) and some of these may not be 
solution to the GF equations. Our numerical experience 
indicates that this is not a practical problem (ie we always 
find a GF solution), but we note the technical possibility 
for completeness. Further, in the next section, we resolve 
this by picking particular choices for IV that eliminate this 
possibility. 




Theorem III. lisa general result that describes a family of 
monotonicity domains, parameterized by W. However, there 
are two drawbacks of this approach: 


• It is unclear how to choose W. Ideally, we would like 
to choose W so that C (W) is as large as possible and 
includes all the practically relevant solutions (solutions 
satisfying typical bounds on flows that one expects will 
hold in practice). 

• The constraint ([S) does not have a simple interpretation, 
so it is difficult to relate the constraints to typical 
operational constraints imposed on gas flows. 

In Section [III-A| we describe how to choose IT so that the 
LMI based condition (| 8 ]l simplifies to a condition involving 
only bounds on the flows. 


A. Selection of Monotonicity Domain 

In this Section, we show that for particular choices of 
IT, the LMI condition ([8]) simplifies significantly. Theorem 
|III.2| shows that if a condition relating the matrices A a 
and A is satisfied, then the GF operator is monotone for 
ip > 0. The basic intuition is that since matrix in |8]i has 
constant terms A , M, these can be cancelled out to get 0 by 
appropriate choice of IT. Similarly, one can cancel out the 
terms involving <p as well, so the overall condition reduces 
to a simple criterion on ip. 

Our first result shows that for the special case of trees and 
systems with no compression, which are already known to 
be efficiently solvable [?][?], the GF operator is monotone 
as long as ip > 0. Thus, any GF solution must lie within 
the monotonicity domain of the GF operator and can be 
efficiently found. Thus, our approach recovers the previously 
known results as special cases. 

Theorem III.2 (Exactness for Trees and Systems with No 
Compression). Suppose 

Aj (a t (AA t ) - 1 A - /) = 0 (9) 

This is guaranteed if one of the following conditions is 
satisfied: 

• The network is a tree. 

• There are no compressors, that is, ak = 1 for all k £ £ 
Then, with the choice 

(AjA T (AA T y 1 0 0 \ 

W = I 0 di (b) 0 

V 0 0 di (b) j 

WF (Cf) is a monotone operator over the domain ip > 0. Let 
pi > 0 be any positive number. Define the domain: 

Cp = {C:^>0,^>p,A^ = q} 

Let (f be any solution to the monotone VI: 

c G Cp,{WF(Q,x-0 > 0 VxGCp 

If |Q = Q. cmd there exists 7r > 0 satisfying 

A a n = b*Q*Q 


then yir, Q, Cfj is a solution to the GF equations. If not, 
there are no solutions to the GF equations. 

([9]) is a very special condition that is only satisfied for 
specific choice of a for networks with general topologies. 
Hence the above theorem has limited applicability for general 
networks. The following theorem gives a way of choosing 
IT such that any solution to the GF with certain upper and 
lower bounds on ip can be found efficiently. 

Theorem III.3 (Flow bounds for general system). Let 
W a £ R raxn ,fr 6 £ R mxm ,7 be an optimal solution to 
the following quasi-convex optimization problem 


Maximize 7 
Subject to 

( di (jp) W b - di (IT C )\ 

( v (fr b -di(ir c )) T ^di(ir c ) ) y0 

W a A = (w b di{bf 1 A a y 


Sy (W b )..-*> 7 (W b ) 

\j¥* 

Then, with the choice 


33 


GO) 

( 11 ) 

( 12 ) 

(13) 

(14) 


/IT° 0 0 

W= 0 W b 0 

\ 0 0 di {W c )j 

For any (3 > 0, WF iff) is a monotone operator over the 
domain C ( y 7 given by: 


{C : C* > 0,141 < C/M < 7/3} 


Let be any solution to the monotone VI: 


t£Cp n ,(WF(0,x-() >0 Vx£C^ (15) 
if ic;i = c; and there exists 7r > 0 satisfying 
Aa 7T = b*Q*Q 

then ^ 7 T, QiCj) ls rf ie unique solution to the GF equations 
in C/ 3 :1 . If not, there are no solutions to the GF equations in 

C/3,7- 


IV. Numerical Examples 

We test our results on the simplest non-tree network: A 
Kite network shown in Fig. [T] 


4 



Fig. 1: Kite Network 








For this network, we find a value of 7 = 54. We observe 
that for this injection configuration, the flows are well within 
the bounds imposed by 7 . 



Fig. 2: Synthetic 16 bus network 


We also tested our approach on a larger 16 bus network 
shown in Fig. [2] We consider three configurations of com¬ 
pressors in the network generated for three different patterns 
of injection and consumption at the sources and sinks respec¬ 
tively. For these three configurations of compression, we find 
the values of 7 to be 470,185,85 respectively. We observe 
that all the numbers are well above the actual bounds on the 
GF. Thus, in all these cases, the GF equations can be solved 
by solving a monotone variational inequality ( fT5) >. 

Numerically, it seems like the value of 7 reduces as the 
mean compression in the network increases. To study this 
effect, we scaled up the compression ratios uniformly and 
studied the change in 7 . The results are illustrated in Fig. 0 



Fig. 3: Effect of Compression Ratios on 7 


equations lie at the heart of many computations in the natural 
gas network studies related to operations, control, optimiza¬ 
tion and planning. Therefore, to solve the GF equations in the 
practical cases of networks, is an important problem which, 
to the best of authors’ knowledge, was not given much of 
attention so far. 

In this manuscript we remedy the problem and present 
a novel approach based on monotone operator theory to 
solve the GF equations. We characterized a whole family of 
monotonicity domains of the GF equations. Within each of 
these domains, determining existence of solutions and finding 
solutions to the GF equations is easy and can be done by 
solving a monotone variational inequality. Further, we show 
that by solving a quasi convex optimization problem, one can 
find a monotonicity domain that contains all GF solutions 
satisfying certain flow bounds. Numerical studies show that 
these bounds are sufficiently tight to capture all solutions 
of practical interest. Thus, for all practical purposes, our 
approach can efficiently find (or prove non-existence of) 
solutions to the GF equations. In future work, we plan to 
study if the monotone operator approach developed in the 
manuscript can be extended to prove that there is always a 
unique solution to the GF equations (this is indeed what is 
observed in our experiments). We plan to focus in the future 
on developing efficient scalable algorithms to solve large GF 
problems using the monotone operator approach. Another 
direction for future analysis is to extend the monotone 
operator approach to more general dynamic case, of interest 
in short term (hours and beyond) operations. Finally, we 
also envision that our approach may be critical for boosting 
performance of more complex multi-level optimization and 
control problems where GF are embedded into low-level 
condition(s). 
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Fig.0 shows that the value of 7 is fairly large (above 10 ) 
for compression ratios of up to about 7, which is already well 
above possible practical value. Hence, we believe that, for 
all practical purposes, our approach can always solve the GF 
equations (or certify non-existence of practical solutions). 

V. Conclusion 

As the interdependency between gas and power systems 
increases, it will become critical to have efficient and ac¬ 
curate tools for gas system operations. The steady state GF 


VI. Appendix 
A. Proof of theorem m 

Proof The Jacobian of the gas flow operator F is given by 

/ 0 A 0 \ 

I A a — di (6 * #) —di (b * <j)) ) = 

\0 di (#>) -di (#) ) 

(I 0 0 \ / 0 

0 — di (6) 0 j -di (b)~ 1 A a 

\o 0 -17 V 0 


A 0 \ 
di (#>) di ( 0 ) 

—di (4>) di (#)/ 
















Let ( £ C (W). It is easy to see that 

0 


1 ° 

Sy ( W 0 -di(6) _1 

\0 0 

Thus, the operator defined by 


0 | VF(C) | FO 


B. Proof of theorem III. 2 


Proof With this choice of W, the monotonicity domain 
becomes condition (|8]i becomes 

0 W a A 0 

Sy I I — A a di (b * ip) di (b * cp) 

0 —di (b * </>) di (b*ip). 


A 0 


The matrix evaluates to 

0 W a A - Aj 

2 di (b * ip) 


W a A r - A a 

0 


-0 


0 
0 

2 di (b * ip)j 


A 0 


If W a A = A C , T . then this condition is always true when ip > 
0. However, W a A = A a T is an overdetermined equation in 
general since W a £ R n x " and the number of equations 
is nm (A a > n for a connected network). However, if a 
solution exists, it must also be a solution to 

W a AA r = A a T A T =>- W a = A a T A T (AA T y 1 

Plugging this back into the original equation, we get 

A a T A T [A A T y 1 A = A a T 

=> Aj (A T (AA T ) _1 A ~l)= 0 

This is exactly the condition of the theorem. 

From the definition of A a , A a = A 7 if a = 1, so that 


A n 


( A t (AA' 1 


, -1 


A- I) = 


aa t (aa 7 


A~ A = A-A = 0 


Now, we can write 

F w (x) - F w (x ) = 


rl dF ( tx + (1 — t)x) 


dr 


dr 


= f VF (tx + (1 — t)x) {x — x) dr 

Jo 


(I 

0 

°\ 


0 

-di(fe)" 1 

0 F (a:) 

= [ 

Vo 

0 

-I 

Jo 


satisfies Sy (VG (x)) A 0. Thus G ( x) is monotone over the 
domain x £ C {W) and C, (W) is a monotonicity domain of 
F. □ 


Since the bottom 2m x 2 to block of 
Sy (VF (tx + (1 — t)x)) is positive definite, the integral 
is non-zero unless x^ = x^^x^ = x Thus, if there is a 
solution to the GF equations in Cp (which is also a solution 
to the VI), then given any solution Cf of the VI, we must 
have IQI = c; and there must exist n > 0 satisfying 

A a w = b*Q*Q 

Otherwise, there are no solutions to the GF equations in 

C/3. □ 


C. Proof of theorem III.3 


Proof. For this choice W, the condition <[8]> evaluates to 


Sy ( ( -W b di(bT x A 0 


W a A 0 

W b di (ip) W b di (cp) 

— &i(W c *cp) di (W c *ip); 


being positive semidefinite. The above matrix evaluates to 


0 Q 

Q T Sy (W b di (ip)) 

0 di(cP)(W b -di(W c )) J 


0 


(W b - di (W c )) di (cp) 
2di (W c * ip) 


where Q = W a A — A a T di (b)~ W b . From the conditions 
of the theorem, we have Q = 0. Hence, we only need to 
prove that 

Sy (W b di (iP)) (W b - di (W c )) di (f)\ ^ 
di (tp) (W b — di (W C )) T 2di (W c * ip) )~ 

Since ip > 0, this condition is equivalent to (using Schur 
complements) 


If the network is a tree, A is an invertible square matrix, 
so that 

A T (AA T y 1 A = A T (A T y 1 (A)~ 1 A = I 

Hence the condition is satisfied in this case as well. 

Uniqueness of solution to the VI: The key observation 
is that for ip £ Cp, F\y is strictly monotone in its last 
2 to coordinates, since the bottom 2m x 2m sub-block of 
Sy (VFvk) is positive definite. It is easy to see that if there 
are two solutions a;, a; to the VI, they must satisfy: 

(Fw (x) — Fw (x),x — x) = 0 


2Sy (W b d\ (ip)) A 
(W b - di (W c )) di 


(W b -di(W c )Y 


^ W c * ip y 

Since \<p\ < ip, tp 2 < ip, so the above condition becomes 
2Sy (W b d\ (ip)) 

A (W b - di (W c )) di (W b - di (W C )) T 

This is true if 

2Sy ( W b d\(iP)) 

A (W b - di (W c )) di (W b - di (W C )) T 






A sufficient condition for the above matrix to be positive 
definite is 

2di (77) /? y 

(W b - di (W c )) di ( w b - di ( W C )) T 

(Sy(^ 6 ))«-^>7^El(Sy(^ 6 ))«lj 

which is equivalent to 

( 2di (77) (fU b -di(fU c ))\ 

l y (W /6 -di(W /c )) T ±di (W c ) ) ~ 

(Sy(^)) ii -77 i > 7 ^|( s y(M/ b )) ij |j 

This is exactly the condition that is required of W in the 
theorem. 

Uniqueness follows from a similar argument as the theo¬ 
rem above. □ 



